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(Dated: September 3, 2010) 

In 1998, Shibata and Sasaki [Phys. Rev. D 58, 104011 (1998)] presented an approximate analytical 
formula for the radius of the innermost stable circular orbit (ISCO) of a neutral test particle around 
a massive, rotating and deformed source. In the present paper, we generalize their expression by 
including the magnetic dipole moment. We show that our approximate analytical formulas are 
accurate enough by comparing them with the six-parametric exact solution calculated by Pachon 
et. al. [Phys. Rev. D 73, 104038 (2006)] along with the numerical data presented by Berti and 
Stergioulas [MNRAS 350, 1416 (2004)] for realistic neutron stars. As a main result, we find that in 
general, the radius at ISCO exhibits a decreasing behavior with increasing magnetic field. However, 
for magnetic fields below 100GT the variation of the radius at ISCO is negligible and hence the 
non-magnetized approximate expression can be used. In addition, we derive approximate analytical 
formulas for angular velocity, energy and angular momentum of the test particle at ISCO. 

PACS numbers: 04.40.Nr, 04.25.Nx, 04.40.Dg, 04.20.Jb 



I. INTRODUCTION 



The discovery of quasi-periodic oscillations (QPOs) with 
frequencies around 1 kHz from several low-mass X-ray 
binaries (LMXBs) [1] has been increasing the interest in 
the detailed theory of disk accretion onto neutron stars. 
Several authors have suggested that at least some of the 
kHz QPOs may be related to the Kepler frequency at the 
innermost stable circular orbit (ISCO) of the accretion 
disk around a neutron star (see e.g. [1, 2]). Stergioulas 
et al. [3] have suggested that the frequency of the co- 
rotating orbit at ISCO in compact stellar remnant could 
be determined through X-ray observations of low-mass 
X-ray binaries and it could be used to constraint the 
equation of state (EOS) of ultradense matter. Morsink 
and Stella [4] have remarked the central role of ISCO in 
relativistic precession of orbits around neutron stars and 
Bulik et al. [5] have shown that observations are consis- 
tent with the assumption that the maximum-frequency 
QPOs occurs at the ISCO. The last statement could be 
used to test general relativity (GR) in the strong-field 
regime around accreting neutron stars, or even to mea- 
sure the stellar mass by directly comparing the highest 
frequency manifest in the X-ray flux with the relativistic 
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formula for the orbital frequency in the ISCO orbit [6]. 

On the other hand, the study of the structure and dynam- 
ics of neutron stars endowed with magnetic field in GR 
is an active, interesting and challenging theoretical issue. 
The influence of magnetic field on the properties of neu- 
tron star rotating at the Kepler frequency has been shown 
in Ref. [7] . In Ref. [8] , Broderick et al. have studied the 
implications of very strong magnetic fields on the struc- 
ture of neutron stars; in particular, Cardall et. al. [9] 
have indicated how magnetic field affects the maximum 
mass of stars. In Ref. [10], the ellipticity of the deformed 
star due to the rotation and to the magnetic field is cal- 
culated, and these two effects arc compared to each other 
within GR. In addition, the formulation of deformation 
of relativistic stars due to the magnetic stress, consid- 
ering the magnetic fields as perturbations from spherical 
stars, has also been studied in [11] by means of an analyt- 
ical treatment assuming weak magnetic fields compared 
to gravity. The quadrupole deformation of magnetized 
Newtonian stars was discussed by Chandrasckhar and 
Fermi [12] and Ferraro [13]. The GR approach was done 
fully numerical by Bonazzola and Gourgoulhon [14] and 
Bocquct et al. [15], who pointed out that deformations of 
the stars induced by magnetic fields become appreciable 
only for fields greater than 10GT. 

More than one decade ago, Shibata and Sasaki [16] (here- 
after S&S) computed an approximate analytical formula 
for the radius at ISCO on massive rotating and arbitrarily 
deformed sources within GR. They considered the role of 
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the quadrupolc moment of mass in physics related to neu- 
tron stars (this fact has also been noted by other authors, 
see for instance [17] and references therein), by including 
multipolar moments of mass up to the 2 4 -pole order in 
their calculations. Yet, as discussed above, there exist an 
strong influence of the magnetic field on the structure of 
the neutron stars and hence is desirable to include also 
the magnetic field in an approximation as the given by 
S&S. In the present paper we tackle this point by follow- 
ing the procedure by S&S to compute approximate formu- 
las for the ISCO including the magnetic dipole moment. 
Thereby, our goal is to calculate approximate formulas for 
the radius, angular velocity, energy and angular momen- 
tum at ISCO for massive rotating and deformed sources 
endowed with magnetic dipole. Due to the importance 
that these parameters have for magnetized neutron stars 
(see e.g. [15, 18]), we should take into account in our 
treatment at least the physical parameters of mass, an- 
gular momentum, quadrupolar moment of mass, current 
octupole moment and magnetic dipole. 

The plan of this paper is as follows: In Section II the gen- 
eral formalism to calculate the ISCO for a neutral parti- 
cle orbiting around a massive source in GR is presented. 
The procedure to compute the approximate formulas for 
radius , angular velocity, energy and angular momentum 
at ISCO of a neutral test particle is shown in Section 
III. The results along with their analysis are presented 
in Section IV. Finally, we present the conclusions of our 
study. 



II. ISCO AND THE MULTIPOLAR 
STRUCTURE 



The metric describing the geometry of space-time around 
a stationary and axisymmetric source, can be written as 
[191 



ds 2 



f(dt- ud4) 2 + f ~ 1 [e 27 {dp 2 + dz 2 ) + p 1 d<f> 2 } , ( 1 ) 



where /, 7 and u are functions of the quasi-cylindrical 
Weyl-Lewis-Papapetrou coordinates (t,p,z,(f)). In this 
paper we use geometrized units c = G = 1 . Hence all the 
physical quantities are measured in units of length [L] . 

In a standard way, we use the line element (1) to find 
the geodesic equations for a neutral test particle on the 
equatorial plane, which reads as follows: 



rfi _ Eg vv + Lg tip 
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with g tt= = -/, g tv = fu and g vv = -fuj 2 + p 2 /f, and 
V e ff(p) denotes the effective potential. 

Circular prograde (or co-rotating) orbits will occur at 
radius p when V e // = and dV e ff/dp = 0, which im- 
poses the following conditions for the angular velocity f2, 
the energy E and the angular momentum L of the test 
particle, 
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The stability of the circular orbit is determined by the 
sign of 
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hereby, ISCOs will occur if and only if d 2 V e ff/dp 2 = 0. 
It is worth mentioning that the above formulae do not 
depend on the metric function 7 and therefore we will 
left aside this metric function in the rest of the paper. 

In order to calculate the metric functions in the elec- 
trovacuum case, we use the Ernst formulation [20]. Via 
Ernst's procedure, the Einstein-Maxwell equations can 
be reformulated in terms of the complex potentials £ and 
$ as 



(Re(£) + |$| 2 )V 2 £ 
(Re(£) + |$| 2 )V 2 $ 



(V£ + 2<TV$) • V£, 

(V£ + 2$*V$) • V$. (8) 



Once the potentials are known, the metric functions / 
and us can be constructed by using 



£ = 



/-|$| 2 +ilm(£), 
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dz 
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dz 
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(10) 



For getting a more intuitive and physical approach, is 
helpful to change the potentials £ and $ to the potentials 
£ and q throughout the following definitions 



£ 



i + e 



(ii) 



This change elucidates the procedure because the poten- 
tials £ and q are related to the gravitational and electro- 
magnetic moment of the source in a very direct way. In 
order to calculate the multipolar moments of an asymp- 
totically flat space-time, according to the Geroch-Hanscn 
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definition [21, 22] we use the procedure of Fodor et al. 
[23] with the corrections given by Sotoriou and Aposto- 
latos [24]. We need to map the initial 3-metric to a con- 
formal one hij hij = ft 2 hij. The conformal factor f2 
should satisfy the following conditions: Q|a = DiCl\\ = 
and L>iDj£l\\ = 2/i^ Ia, where A is the point added to the 
initial manifold that represents infinity. f2 transforms the 
complex gravitational and electromagnetic potentials £ 
and q into £ = il _1 / 2 £ and q = Q~ x / 2 q respectively. The 
conformal factor is given by SI — f 2 — p 2 + z 2 , and the 
transformation between unbarred and barred variables 
reads as 

P=—> 2> E= ~, 2' = & ( 12 ) 

p l + z 2 - p l + Z A 

I 



which brings infinity at the origin of the axes (p, z) = 
(0, 0). The potentials £ and q can be written in a power 
series expansion of p and z as 



i=Y f a ij p i z^, q=Y i b ij (?zt. (13) 

i.J—O i-J—0 

Due to the analyticity of the potentials at the axis of 
symmetry. Oy and bij must vanish when i is odd. The 
coefficients in the above power series can be calculated 
by using the relations [24] 



(r + 2) 2 a r+2 . s = -(s + 2)(s + l)a r , s+2 + ^ {a H a* mn - b H b* mn )[a pg (p 2 + g 2 - 4p - 5g - 2pk - 2gl - 2) 
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and 
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+ b p+ 2, g -2(p + 2)(p + 2-2k) + b p - 2 , g+ 2(g + 2)(g + l-2l)), (15) 



where m — r — fc — p, < k < r, < p < r — k with 
k and p even, and n = s — I — g, < / < s + 1, with 
— 1 < g < s — I. These recurrence relations could build 
the whole power series of £ and q from their values on 
the axis of symmetry 

00 00 
£(p = 0,z) = J2 m iZ i > ?(P = 0,f) = X)«« i . ( 16 ) 

i=0 i=0 

where the coefficients in the above series expansion are 
related to the values of the multipolc moments of the 
space-time g, = boi and = aoi [23, 24]. 



III. APPROXIMATE FORMULAS AT ISCO 



Following the scheme given by S&S [16] we assume that 
the only non vanishing multipolc moments are the mass 
Mo = M, the angular momentum M\ = qM 2 , the mass 
quadrupole A/ 2 = — Q2M 3 , the current octupole mo- 
ment M 3 = -Q3M 4 , the 2 4 -pole M 4 = Q 4 M 5 and 
additionally, the magnetic dipole moment Ai = pM 2 , 
where q, Q2, Q3, Q4 and p are dimcnsionlcss parameters. 
In order to keep the approximation consistent up to 



C(e 4 ), we formally set q -> eq, Q 2 -> e 2 Q 2 , Q3 -t 
e 3 Q3, Qi £ 4 Q± and p -t e 2 p. 

With the aim to calculate the approximate potentials (13) 
as a truncated power series, we carry out the following 
steps: (i) We compute the gravitational and electromag- 
netic multipolcs up to order twelve (using the corrected 
formulas given by Sotiriou and Apostolatos [24]), as a 
function of the coefficients on the symmetry axis, a j 
and 60 j . (ii) By inverting these expressions, we then get 
the coefficients oqj and boj as a function of the multi- 
poles (see the Appendix and the note [25] at the end of 
the paper), (iii) Then, we use the expressions for aoj 
and bo.j along with the recurrence relations (14)- (15), in 
order to calculate the coefficients aij and bij up to 0(e 4 ). 
We do not present here these quantities because of their 
cumbersome form, but they are available under request 
to the authors. 

Once we know the approximate expressions for the com- 
plex potentials (13), it is possible to obtain the approxi- 
mate expressions of the Ernst potentials £ and $ by ap- 
plying Eq. (11). Consequently, we compute the metric 
functions / and uj (9)-(10), by expanding in power series 
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of the inverse of p: 
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fe=i 



k=l 



(17) 



(18) 



-^ISCO 

M 



3.46410 - 0.94281? - 0.44452/ + 0.18793Q 2 

- 0.39579/ + 0.29982gQ 2 

- 0.00519g/ - 0.44854/ 
-0.05055Q 2 . - 0.09282 g Q 3 + 0.00935Q 4 

- 0.49506g 2 Q 2 - 0.01099g 2 / 
(23) 



-0.00195/ - 
-0.03926Q 3 



-0.00170Q 2 ^ 
-0.00002^. 



where Cf >k and C^^k are functions of the multipolcs. 

By using the Eqs. (17) and (18) we can cast the Eq. (7) 
as 



IV. RESULTS AND ANALYSIS 



e k A k ( Pl M, q, Q 2 , fis, 24, A*) = 0. (19) 



k=0 



Solving it for p, we obtain for the circumferencial radius 
R = at ISCO: 
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1 - 0.54433g - 0.22651/ + 0.17992Q 2 

-0.00323/ - 0.23122/ + 0.26353gQ 2 
-0.05318Q3 - 0.00765?/ - 0.29981/ 
+0.44887/Q 2 - 0.06260Q 2 - 0.11325gQ 3 
+0.01546Q 4 - 0.01572// + 0.00312Q 2 / 



-0.00004/ 



(20) 



and for the angular velocity, the energy and the angular 
momentum at ISCO: 
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1 



6V6M 
+0.00446/ 
+0.07433 Q 3 
+0.12813Q 2 . 



(1 + 0.74860q + 0.78106/ - 0.23432Q 2 
0.98328/ - 0.64492gQ 2 



0.01608g/ - 
0.25050gQ 3 



1.38828/ 
-0.02132Q4 



-Eisco 



-0.00596Q 2 / - 1.4235l/Q 2 + 0.04191// 
+0.00007/), (21) 

0.94280 - 0.03208q - 0.02977/ + 0.00794Q 2 
-0.00010/ - 0.03417/ + 0.01980gQ 2 
-0.00200Q 3 - 0.00035q/ - 0.04427/ 

0.00049Q 4 



-0.00331Q 2 . - 0.00621qQ3 



-O.OOOI2Q2/ 
-O.OOOOl/, 



0.04044/ Q 2 - 0. 



V 

(22) 



First, we verify our results in the vacuum case. We 
check our formula for the radius at ISCO (20), by set- 
ting Q 2 = /, Q 3 = /, Qi = q A and p = 0, and compar- 
ing it with the exact expression given by Bardeen et al. 
[26]. Here, we find that the error for q < 0.5 is smaller 
than 1%. Next, we compare our expression (20) with the 
six-parametric exact solution given by Pachon et al. [27] 
along with the numerical data given by Berti and Stcr- 
gioulas (hereafter B&S) for selected EOS [28]. In this 
case, for q < 0.3, the error still being smaller than 0.6% 
in all the cases. It should be noted at this stage that de- 
spite the mistake in Eqs. (2.24), (2.25), (Bl) and (B2) of 
S&S (see note [25]), the difference with our expressions 
(20)-(23) is smaller than 0.1% and therefore negligible. 

Now, we turn our attention to the electrovacuum case. 
Starting from the following limiting statements: (i) Mag- 
netars lose rotational speed very quickly due to their 
high magnetic field, (ii) Given their rarity, the possi- 
bility to observe a new-born rapidly rotating magnctar 
is negligible, (iii) The amount of data for observed mag- 
netars is minimum (cf. [29]) without any data of their 
higher-order multipolar structure, (iv) Theoretical stud- 
ies reporting numerical data of the multipolar structure 
of magnetars (cf. [15]), did not present numerical data 
for the radius at ISCO nor numerical data for higher 
mass-rotation multipole moments, (v) The Pachon et 
al. solution [27] fits very well with realistic numerical 
interior solutions for slowly rotating neutron stars, pos- 
sessing an arbitrary magnetic dipole parameter. Let us 
to assume that the Pachon et al. solution is a good model 
for realistic slowly rotating magnetars. 

With the aim of testing our approximate formulae, we use 
the parameters calculated by B&S for neutron stars, ex- 
trapolated to the case of a non- vanishing magnetic dipole. 
In order to observe the effect of the magnetic field on the 
radius, hereafter we restrict ourselves to the use of the ap- 
proximate formula (20) in the presence of magnetic dipole 
p 7^ 0, where this approach will be labeled as "Mag" or 
in its absence p = 0, where this approach will be labeled 
as "Non-Mag". 
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A rough estimate of the numerical solutions of the 
Einstein-Maxwell equations presented by Bocquet et al. 
(see Table 2 in [15]) for models of rapidly rotating mag- 
netized neutron stars, suggest that the magnetic dipole 
M. belongs, in average, to the interval to 10 32 Am 2 , 
corresponding to magnetic fields in the range to 10 12 T. 
Hence, from the proportionality between M. and the mag- 
netic field B, the observed magnetic field for magnetars 
ca. 10 11 T [29], should roughly correspond to a magnetic 
dipole moment of about 10 31 Am 2 . Moreover, in natural 
units the magnetic dipole moment has an order of 



(a) 



(b) 



M = Y- — Ms.i, 



(24) 



which represents a value of the magnetic dipole of Ai ~ 
lkm 2 for observed magnetars (see for instance the para- 
graph around Eq. (20) of Ref . [27]). 

For the particular case M. = lkm 2 , in Fig.l we plot 
the radii at IS CO for the six parametric exact solution 
presented by Pachon et al. (Six-Parametric), the ap- 
proximate formula (20) with magnetic dipole (Mag) and 
the approximate formula without magnetic dipole (Non- 
Mag). Here, we use the data of B&S for EOS AU se- 
quence with constant rest mass corresponding to a non- 
rotating model of 1.578M (panel a) and the EOS APRb 
sequence with constant rest mass corresponding to the 
maximum- mass model in the non-rotating limit 2.672M Q 
(panel b). As can be easily noted from the insets in Fig.l, 
for the case of a dipolar magnetic moment of 10 31 Am 2 , 
the changes introduced in the radius at ISCO by the ap- 
proximate formula (20) taking into account the magnetic 
dipole are negligible. 

In Fig. 2 we plot the radii at ISCO for the same param- 
eters as in Fig.l, but using the high value M. = 10km 2 . 
For EOS AU (panel a) the error for the case q max ~ 0.7 is 
around 4% for the Non-Mag. case, while it is just close 
to 1% for Mag. For EOS APRb (panel b) the errors are 
1% for Non-Mag. and 0.5% for Mag. Therefore, it can 
be concluded that for very strong magnetic fields (ten 
times larger than the ones observed for magnetars) the 
influence of the magnetic field on the radius at ISCO is 
significantly important, producing a decreasing tendency 
on the radius at ISCO with the increasing of the magnetic 
dipole strength. 

Finally, just to illustrate the influence of the magnetic 
dipole M. on the radius at ISCO for magnetars, in Fig. 3 
we plot the radii at ISCO for the six-parametric exact 
solution, the magnetized approximate solution and the 
non-magnetized one. Here we set the multipole moments 
in accordance with typical values for observed magne- 
tars: M = 1.402M Q , q = Q 2 = Q 3 = Q± = 10 -3 . It 
can be seen that the radius at ISCO is affected by the 
existence of magnetic dipole, with an appreciable varia- 
tion with respect to the non-magnetized case. At \i ~ 2 
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FIG. 1: Case M aa t = lkm 2 : Radius at ISCO (henceforth 
measured in Km) for the EOS AU sequence with constant 
rest mass corresponding to a non-rotating model of 1.578-Mq 
for prograde orbits (panel a) and for EOS APRb sequence 
with constant rest mass corresponding to a sequence that ter- 
minates at the maximum-mass non-rotating model in the 
non-rotating limit of 2.672A/0 for prograde orbits (panel b). 
The difference between the Mag. and Non-Mag. cases is neg- 
ligible, as depicted in the enlargements (insets). Color insets 
online. 



(a) 



(b) 



13 



10 



1 1 r 

Six-Parametric — 

Non Mag. 

Mag. — 




J l l_ 



0.2 0.4 0.6 O.i 
j 




0.2 0.4 0.6 0.8 
j 



FIG. 2: Case M = 10km 2 : Same EOS as in Fig.l. The differ- 
ence between the Mag. and Non-Mag. cases is appreciably 
more marked. Color insets online. 



(corresponding to a value of the magnetic dipole of about 
0.5 x 10 32 Am 2 ), the errors are close to 3% for the non- 
magnetized version and around 1% for the magnetized 
version. With this example we intent to show that the 
same tendency as discussed above holds for realistic val- 
ues of magnetars. 

From a mathematical point of view, the strong influence 
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eluding the magnetic dipole moment of the source), the 
radius at ISCO exhibits appreciable deviations with a ten- 
dency to decrease away from its non-magnetized version, 
only if the magnetic field strength is higher than 100GT. 
Finally, we want to remark that theoretical models for 
slowly rotating magnetars including numerical data for 
the higher multipolar structure and the radius at ISCO 
are desirable. 
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FIG. 3: Dependence of i?isco with fj, for parameters of ob- 
served magnetars. Note that when /i becomes larger the ra- 
dius becomes smaller and the formula with magnetic dipole 
(labeled as Mag.) is close to the analytical solution (labeled 
as Six-parametric.) 



of the magnetic dipole for fi > 1, can be explained tak- 
ing into account that the magnetic dipole parameter /i 
appears into the metric functions by means of |$ | and 
products of imaginary and real parts of $, so for the di- 
mensionless magnetic dipole parameter \x its lower order 
in the metric functions and therefore in the ISCO formu- 
las is of the order /i 2 . From a physical point of view, it 
is a well known fact that in general relativity the mag- 
netic field could induce observable effects into the space- 
time (see e.g. [30] and references therein) and therefore 
changes onto the dynamics of non magnetized neutral 
particles orbiting around such sources can be expected 
[311. 



CONCLUDING REMARKS 



In this paper, we have obtained simple approximate for- 
mulas for the radius, the angular velocity, the energy and 
the angular momentum of a neutral test particle at the 
innermost stable circular orbit (ISCO). These formulas 
correct the results obtained by Shibata and Sasaki [16] 
and generalize it to the case of neutral test particles mov- 
ing on the equatorial plane around a rotating source en- 
dowed with magnetic dipole. In order to test the accuracy 
of our approximate expressions, we first have compared it 
with the radius at ISCO (calculated for non-magnetized 
neutron stars) by means of the Kerr solution, the six- 
parametric exact solution given by Pachon et al. [27] and 
the numerical data given by Berti and Stergioulas [28] for 
selected EOS. In all cases, our formula differs from the 
numerical results by at most 0.6% in the slow-rotation 
regime (i.e. for j < 0.3). As a main result, it was found 
that when using realistic parameters for magnetars (in- 
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APPENDIX: COEFFICIENTS FOR THE ERNST 
POTENTIALS ON THE SYMMETRY AXIS. 



The following expressions show the data on the symmetry 
axis mi := ao.i and qi := bo.i, obtained by using the 
procedure of FHP [23] with the formulas given by Sotiriou 
and Apostolatos [24], in order to get a source with the 
following multipole structure: Mass Mq = M, angular 
momentum Mi = qM 2 , mass quadrupole M 2 = — Q 2 M 3 , 
current octupole moment M3 = — Q3M 4 , 2 4 -pole of mass 
AI4 = Q4M 5 and magnetic dipole moment A4 = ixM 2 
(all the other multipolar moments are set to zero), 



a 0i0 = M, a A = iM 2 q, 



ao,2 = ~M Q 2 , a , 3 
1 



-iM 4 Q 3 , 



a 0A = -^M 5 (lOq 2 - 3/j 2 - 10Q 2 + 70Q 4 ) , 



a ,5 = -i^M 6 (q 3 - 1^1 - 8Q 2 q + 7Q 3 ) , 



ao,6 - M 



^ 231 21 33 77 
205Q 2 /i 2 32 M 2 Q 2 , 6Q 4 \ 



4158 



945 21 11 
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ao,7 



= iM 8 - 



13g 3 46n 2 q \QQ 2 q 



231 

+Af 8 0(e 5 ), 



M 9 - 



10395 



77 



5Q 3 
33 



3003 



2623mV 536Q 2 g 2 5g 2 
^ 231 



216216 



3003 



140Q 3 g 257^ 4 461Q 2 137Q 2 /i 2 



429 
985^ 2 



216216 3003 
5Q 2 , 45Q 4 



36036 231 



143 



90090 
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(A.25) 
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